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Research  Objectives: 

The  objective  of  our  original  proposal  was  to  use  the  Minimuin  Current  Corona  (MCC) 
model  to  determine  where  strong  currents  are  likely  to  develop  above  an  active  region,  and 
compare  the  locations  with  observed  flare  signatures.  If  there  is  a  good  correspondence,  we 
can  also  estimate  the  energy  available  for  the  flare  from  the  current  sheet  in  question.  In  the 
MCC  model,  the  magnetic  held  in  the  corona  is  represented  by  the  field  due  to  a  collection 
of  point  sources  on  the  boundary:  a  Magnetic  Charge  Topology  (MCT)  model.  Evolving 
these  sources  consistent  with  the  evolution  of  the  boundary  causes  currents  to  begin  flowing 
in  the  corona.  The  minimum  energy  state  occurs  when  currents  flow  only  along  separator  field 
lines:  topologically  significant  field  lines  which  connect  magnetic  null  points.  In  the  absence  of 
reconnection,  the  magnitudes  of  these  currents  can  be  estimated,  along  with  the  energy  stored 
by  each. 

The  objectives  of  the  research  were  broadened  to  include  an  analysis  of  the  injection  of 
magnetic  helicity  into  the  solar  corona.  In  brief,  reconnection  is  likely  to  play  a  crucial  role 
in  the  initiation  or  release  of  energy  in  solar  flares,  while  helicity  may  provide  insight  into  the 
mechanism  responsible  for  Coronal  Mass  Ejections.  W’e  have  used  the  partitioning  needed  to 
implement  the  MCC  model  to  understand  the  way  in  which  helicity  is  injected  into  the  solar 
corona. 

Our  effort  was  also  extended  to  justify  the  use  of  the  MCC  model.  In  particular,  the  scien¬ 
tific  community  has  expressed  skepticism  that  the  MCT  potential  field  is  a  good  representation 
of  the  actual  potential  field,  and  a  recent  letter  by  Georgoulis  &  Rust  (2007)  suggested  that  a 
good  parameter  for  forecasting  major  flares  could  be  derived  from  simply  the  MCT  connectivity 
at  a  given  time.  Thus  we  have  spent  time  addressing  both  of  these  concerns. 

Work  Carried  Out: 

Instead  of  solving  the  full-blown  magnetohydrodyiiamic  (MHD)  equations,  the  MCC  models 
the  coronal  magnetic  field  with  a  collection  of  point  sources  on  the  boundary.  In  this  case, 
the  topology  becomes  particularly  simple:  with  a  few  special  exceptions,  field  lines  begin 
and  end  on  sources.  One  of  the  important  exceptions  is  separator  field  lines,  which  begin 
and  end  on  magnetic  null  points,  and  encircle  bundles  of  flux  connecting  pairs  of  sources. 
The  locations  and  magnitudes  of  the  points  sources  (magnetic  charges)  are  determined  by 
partitioning  the  boundary  field  into  flux  concentrations;  each  partition  is  represented  by  a 
single  point  source.  The  potential  field  due  to  the  collection  of  point  sources,  and  thus  the 
amount  of  flux  connecting  each  pair  of  sources,  is  calculated  at  two  times.  The  changes  in  the 
potential  field  connectivity  are  canceled  by  introducing  currents  along  separators,  and  the  free 
energy  introduced  is  computed  from  these  separator  currents.  The  advantages  to  the  model  are 
that  it  shows  where  the  energy  is  stored  and  how  much  can  be  liberated  by  a  single  reconnection 
event,  and  that  it  calculates  the  free  energy  directly,  rather  than  taking  the  difference  between 
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the  total  energy  and  the  potential  energy.  A  disadvantage  is  that  the  model  must  start  from 
a  potential  field  eonfiguration,  and  thus  ean  only  traek  the  evolution  of  the  energy. 

We  have  applied  the  MCC  model  to  multiple  time  series  of  veetor  inagnetograms  from  the 
Imaging  Veetor  Magnetograph  (IVM).  The  evolution  of  the  free  magnetie  energy  was  traeked, 
and  the  loeation  of  the  greatest  energy  storage  was  eompared  to  the  loeation  of  brightening 
in  MCCD  images.  In  addition,  the  MCC  model  was  applied  to  the  mueh  longer  time  series 
of  line  of  sight  magnetograms  available  from  the  Miehelson  Doppler  Imager  (MDI)  on  board 
the  SOHO  spaee  eraft.  The  use  of  these  longer  time  series  required  the  development  of  a  new 
partitioning  algorithm.  The  validity  of  the  partitioning  was  eonfirmed  by  eomparing  the  souree 
veloeities  to  veloeities  determined  using  Loeal  Correlation  Traeking. 

The  new  partitioning  algorithm  allowed  us  to  deeompose  the  magnetie  helieity  injeeted  into 
the  eoroiia  into  a  spin  eomponent  and  a  braiding  eomponent.  The  former  is  likely  to  arise  from 
the  twist  of  sub-photospherie  flux  tubes,  while  the  latter  results  from  motion  of  flux  tubes 
about  one  another.  The  partitioning  algorithm  was  run  for  a  number  of  time  series  of  MDI 
data,  and  used  to  deeompose  the  helieity  injeetion.  The  relative  importanee  of  spin  versus 
braiding  helieity  w^ls  determined. 

Coronal  magnetie  field  ean  be  eharaeterized  by  the  way  its  field  lines  intereoimeet  regions 
of  positive  and  negative  photospherie  flux:  its  eonneetivity.  Conneetivity  ean  be  quantified  as 
the  net  flux  eonneeting  pairs  of  opposing  regions.  It  is  rapid  ehanges  in  this  eonneetivity  that 
are  believed  to  release  stored  magnetie  energy  in  the  form  of  solar  flares.  We  have  performed  a 
detailed  eomparison  of  the  magnetie  field  eonneetivity  for  potential  fields  derived  from  both  an 
MCT  model,  and  from  using  the  observed  magnetie  field  on  the  whole  boundary.  In  addition, 
we  eonsidered  one  further  differenee  in  the  boundary:  the  field  was  either  eontained  within 
eondueting  walls  on  four  sides,  or  was  allowed  to  expand  to  fill  the  half-spaee  above  the 
boundary.  By  doing  this,  we  have  been  able  to  determine  the  limitations  of  determining  the 
eonneetivity  in  an  MCT  model  resulting  from  redueing  the  lower  boundary  to  a  eolleetion  of 
point  sourees,  and  thus  have  established  the  validity  of  this  approximation  in  the  MCC  model. 

Reeent  work  by  Georgoulis  Sz  Rust  (2007)  suggested  that  a  parameter  eonstrueted  from  the 
eonneetivity  of  an  MCT  model  is  “an  effieient  flare-foreeasting  eriterion" .  If  it  were  indeed  the 
ease  that  a  robust  flare  foreeasting  parameter  eould  be  eonstrueted  using  only  the  eonneetivity 
in  an  MCT  model  at  one  time,  then  there  would  be  less  ineentive  to  eonsider  the  MCC  model. 
Thus  we  have  replieated  the  approaeh  of  Georgoulis  Sz  Rust  (2007)  for  our  own  data  base 
of  magnetie  field  observations,  and  eompared  the  performanee  of  their  parameter  to  ones 
previously  investigated. 

Results  Obtained: 

Spatial  Distribution  and  Evolution  of  the  Free  Energy 

We  have  applied  the  MCC  model  to  time  series  of  veetor  magnetograms  from  the  Imaging 
Veetor  Magnetograph  (IVM)  for  several  aetive  regions  whieh  produeed  large  flares.  The  results 
of  the  analysis  for  a  sample  aetive  region  (NOAA  AR 10030)  are  shown  in  Figure  1.  On  the 
left,  the  partitioning  of  the  magnetogram  is  shown,  along  with  the  loeations  of  sourees,  nulls, 
spines  and  fan  traees.  Of  partieular  interest  is  the  eoronal  magnetie  null  at  approximately 
(0,265)  in  the  image. 

To  investigate  the  loeation,  timing,  and  dynamies  of  ehromospherie  and  photospherie  emis¬ 
sion  assoeiated  with  solar  flares,  we  used  data  from  the  Mees  CCD  Imaging  Speetrograph 
(MCCD).  The  MCCD  routinely  obtained  Ha  speetra  with  extremely  high-eadenee  seans  aeross 
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Figure  1:  MCC  analysis  of  NOAA  AR 10030.  Left:  The  vertical  component  of  the  magnetic  field 
(greyscale)  plus  the  partitioning  (lines),  locations  of  sources  (plus  and  cross),  nulls  (triangles,  blue 
for  A- type  and  red  for  B-type,  purple  for  coronal),  spines  (green  lines),  and  fan  traces  (dashed  lines). 
Right:  An  MCCD  image  of  the  thermal  response  in  the  chromosphere  from  Ha  liiie-center  images 
of  the  the  flare  ribbons  which  developed  roughly  10  minutes  after  the  peak  soft  X-ray  flux,  with 
the  same  sources,  nulls,  spines  and  fan  traces.  In  addition,  separator  field  lines  are  shown,  with  the 
color  of  each  separator  determined  by  the  current  along  that  separator,  with  red  being  the  most 
energetic.  Although  the  time  series  used  is  too  short  to  encapsulated  the  full  build  up  of  energy,  it 
is  encouraging  that  some  of  the  most  energetic  separators  are  in  the  vicinity  of  the  brightest  part 
of  the  MCCD  image.  Further,  the  brightening  all  occurs  close  to  or  within  the  separatrix  surface 
of  the  coronal  null  point,  indicating  the  possible  importance  of  such  topological  features.  This  Hare 
occurred  in  conjunction  with  a  Coronal  Mass  Ejection  (CME),  so  the  presence  of  the  coronal  null 
point  is  consistent  with  the  magnetic  breakout  model  for  CME  initiation. 


a  target  active  region.  The  target  for  the  MCCD  was  usually  coordinated  with  that  of  the 
IVM,  making  a  powerful  data  set  of  photospheric  magnetic  field  evolution  and  chromospheric 
thermal  response  and  dynamics  during  solar  energetic  events.  On  the  right  of  Figure  1,  an  Ha 
line  center  image  from  the  MCCD  of  the  flare  is  shown,  with  the  same  sources,  nulls,  spines 
and  fan  traces.  In  addition,  separator  field  lines  are  shown,  with  the  color  of  each  separator 
deteniiined  by  the  current  resulting  from  the  MCC  model  applied  to  a  time  series  of  approx¬ 
imately  one  hour  duration.  Although  the  exact  location  does  not  match,  it  is  encouraging  to 
note  that  several  of  the  most  energetic  separators  are  in  the  vicinity  of  the  brightest  part  of  the 
MCCD  image.  In  fact,  the  chromospheric  brightening  occurs  primarily  within  the  separatrix 
surface  of  the  coronal  null  point.  This  flare  occurred  in  conjunction  with  a  Coronal  Ma^s  Ejec¬ 
tion  (CME),  so  the  presence  of  the  coronal  null  point  is  consistent  with  the  magnetic  breakout 
model  for  CME  initiation. 

The  free  energy  estimate  for  this  active  region  is  much  smaller  than  the  energy  typically 
released  by  a  large  flare.  To  further  irm^stigate  the  build  up  of  energy,  we  have  applied  the 
MCC  model  to  NOAA  AR8210  during  the  time  leading  up  to  the  Ml. 2  flare  on  1998  May  1. 
In  particular,  we  used  both  an  approximately  five  hour  time  series  of  IVM  vector  magnetogram 
data,  and  an  approximately  one  day  time  series  of  MDI  line  of  sight  magnetogram  data.  In 


3 


AW  [erg] 


lime  [h] 


Figure  2:  The  magnetic  free  energy  from  the  MCC  model  for  NOAA  AR8210.  Left:  The  entire 
time  interval  of  MDI  data  considered,  indicating  that  the  energy  builds  up  on  a  timescale  of  a  day  or 
longer.  Right:  Comparison  of  the  energy  derived  from  the  IVM  time  series  (squares)  to  the  MDI  time 
series,  initiated  at  the  same  time  as  the  IVM  time  series  (pluses)  and  initiated  approximately  17  hr 
earlier  (asterisk).  When  initiated  at  the  same  time,  the  resulting  energies  are  very  similar,  indicating 
that  using  the  line  of  sight  component  of  the  field  rather  than  the  vertical  component  does  not  have 
a  large  effect  over  a  short  time  interval,  but  may  diverge  for  longer  times. 


both  cases,  the  currents  along  separators  was  computed,  and  the  free  energy  build  up  inferred. 
Figure  2  shows  the  result  for  the  both  the  full  MDI  time  series  (left)  and  the  shorter  time 
interval  of  the  IVM  data  (right).  It  is  clear  that,  for  this  region,  comparatively  little  energy  is 
injected  during  the  few  hours  immediately  prior  to  the  flare  at  approximately  23  hr.  However, 
over  the  course  of  the  j)revious  day,  a  factor  of  four  more  energy  is  introduced.  Thus  we 
conclude  that,  in  terms  of  the  energetics,  it  is  likely  necessary  to  follow  the  evolution  of  a 
region  for  more  than  a  few  hours. 

There  is  at  present  no  source  of  high  cadence  long  duration  time  series  of  vector  magnetic 
field  obser\^tions  available.  The  long  time  series  are  needed  to  fully  track  the  build  up  of 
magnetic  free  energy,  and  particularly  for  long  time  series,  vector  observations  are  important 
to  separate  out  real  evolution  of  the  solar  magnetic  field  from  evolution  resulting  from  pro¬ 
jections  effects  inherent  to  line  of  sight  observations.  Thus  the  final  results  of  our  work  effort 
were  severely  impacted  by  the  denial  of  our  request  for  a  no  cost  extension  to  our  contract. 
Our  original  proposal  included  applying  the  MCC  model  to  data  from  the  Solar  Dynamics 
Observatory  (SDO),  but  delays  to  launch  of  SDO  mean  they  are  still  unavailable. 

Evolution  of  Magnetic  Properties  for  Two  Active  Regions 

Although  it  is  difficult  to  make  reasonable  energy  estimates  from  the  evolution  of  line  of 
sight  magnetograms,  we  have  applied  the  MCC  model  to  two  active  regions  observed  by  the 
Hinode  X-ray  Telescope  (XRT)  and  the  Transition  Region  and  Coronal  Explorer  (TRACE). 
One  active  region  shows  constant  brightness  in  both  XRT  and  TRACE  observations.  The  other 
active  region  shows  a  brightening  in  the  TRACE  observation  just  after  a  decrease  in  X-ray 
brightness  indicating  the  cooling  of  a  coronal  loop  (see  Fig.  3).  By  considering  the  difference  in 
the  behavior  of  the  two  regions  in  the  MCC  model  (see  Fig.  4),  as  applied  to  MDI  line  of  sight 
magnetograms,  we  hope  to  learn  more  about  changes  in  coronal  energy  release.  This  work  is 
still  proceeding  under  other  sources  of  funding.  An  initial  investigation  based  on  the  difference 
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Figure  3:  The  small  boxes  near  (100"  x  150")  on  top  and  (230"  x  -250")  on  bottom  reprasent  the 
selected  area  of  the  light  curve  for  AR  10938  and  AR  10963,  respectively.  Solid  and  dotted  lines 
represent  the  light  curves  for  XRT  and  TRACE,  respectively.  [Adapted  from  Lee  et  al.  (2009).] 
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Figure  4:  Left:  MDI  magnetogram  with  connectivity  changes  (lines)  and  point  sources  (-!-,x  labeled 
by  numbers).  Thicker  and  darker  lines  show  larger  connectivity  changes.  Right:  Changes  in  the 
Hux  enclosed  by  separators.  Triangles  represent  nulls  -  positive  (red),  negative  (blue),  upright  (moss 
green,  near  180"  x  180"  on  top),  and  coronal  (green,  near  270"  x  —190"  on  bottom).  BlacL  small 
boxes  represent  the  location  for  light  curve  in  Fig.  1.  Yellow  and  orange  lines  represent  larger  Hux 
changes.  The  same  scales  were  used  for  two  active  regions.  [Adapted  from  Lee  et  al.  (2009).] 


between  two  MDI  images  showed  little  difference  between  the  two,  as  would  be  expected  since 
the  time  interval  is  too  short.  Presently  we  are  considering  an  approximately  one  day  time 
interval,  to  see  whether  significant  differences  are  found  in  the  MCC  model. 

Helicity  Decomposition 

Using  the  partitioning  algorithm  developed  for  use  with  long  time  series  of  magnetograms, 
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Figure  5:  Summary  of  the  helicity  flux  (bottom)  and  its  integral  (top)  into  NOAA  AR9114  plotted 
versus  time  in  hours  from  00:00UT  on  Aug-07.  The  total  is  indicated  by  a  solid  line,  its  decomposition 
into  spin  and  braiding  components  is  indicated  by  +  and  x  respectively,  and  the  sum  of  these  is  shown 
as  a  *.  In  this  case,  the  rates  of  injection  of  spin  and  braiding  helicity  are  comparable.  [Adapted  from 
Loiigcope,  Ravindra  &  Barnes  (2007).] 


we  have  developed  a  method  to  usefully  decompose  the  helicity  flux  across  the  photosphere 
into  contributions  of  differing  origin,  called  spin  helicity  and  braiding  helicity.  These  contri¬ 
butions  would  typically  come  at  the  expense  of  twist  and  writhe  helicity,  respectively,  of  a 
sub-photospheric  flux  tube  anchored  to  the  regions.  The  spin  helicity  of  a  given  partition 
quantifies  the  mean  rotation  rate  of  motions  internal  to  that  partition,  while  braiding  helicity 
is  injected  by  the  motions  of  whole  partitions  about  one  another.  This  decomposition  is  shown 
for  NOAA  AR9114  in  Figure  5.  In  this  case,  the  contributions  from  the  spin  and  braiding 
helicity  are  comparable. 

Applying  the  method  to  six  active  regions  shows  cases  where  either  spin  or  braiding  dom¬ 
inates,  and  where  they  have  the  same  signs  and  opposite  signs.  Thus  it  would  seem  that  no 
general  statement  can  be  made  regarding  the  dominance  of  twist  or  writhe  in  supplying  helicity 
to  the  corona. 

Connectivity  Comparison 

To  determine  the  effects  of  the  boundary  conditions  on  the  potential  field  connectivity,  we 
studied  in  detail  one  active  region:  NOAA  AR8636  (Figure  6).  The  active  region  magnetic 
field  was  partitioned  using  the  algorithm  described  in  Barnes,  Longcope  Sz  Leka  (2005),  to 
determine  flux  concentrations.  The  partitioning  was  repeated  for  a  wide  range  of  parameters, 
resulting  in  a  range  in  the  number  of  sources,  from  a  few  dozen  in  the  coarsest  partitioning, 
to  a  few  hundred  in  the  finest.  For  each  way  of  partitioning  the  boundary,  the  connectivity 
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Figure  6:  Connectivity  for  NOAA  AR8636.  The  vertical  magnetic  field  is  shown  on  a  greyscale, 
sources  are  shown  as  -|-,x  for  positive  and  negative  polarity,  and  the  flux  in  each  connection  determines 
the  thickness  and  color  of  the  line  segment  connecting  the  sources,  with  large  flux  connections  shown 
by  dark,  thick  lines.  This  case  has  N  =  116  sources. 


was  calculated  in  four  ways:  using  the  full  magnetogram  as  a  boundary  condition  (M)  versus 
using  point  sources  on  the  lower  boundary  (P),  and  letting  the  field  fill  the  full  half-space 
z  >  0  (H)  versus  bounding  the  field  with  conducting  walls  outside  the  field  of  view  of  the 
magnetogram  (B).  For  one  set  of  partitioning  parameters,  a  comparison  of  the  MB  and  PB 
cases  is  shown  in  Figure  7  (left).  For  large  flux  connections,  meaning  those  with  more  than 
about  10”^  of  the  total  flux  in  the  active  region,  there  is  relatively  little  difference  between  the 
two  cases;  for  small  flux  connections,  the  difference  can  be  large,  with  some  small  connections 
being  completely  absent  in  one  case  or  the  other.  Overall,  there  is  good  agreement  between 
the  connectivities,  with  an  overall  flux  discrepancy  of  6.8%  for  this  case.  That  is,  over  90% 
of  the  flux  in  the  region  connects  the  same  pairs  of  partitions  in  both  ways  of  calculating  the 
potential  field. 

This  result  of  the  large  flux  connections  matching  well,  and  the  small  flux  connections  having 
larger  differences  is  a  general  property  of  the  cases  we  considered.  In  Figure  7  (right),  we  show 
how  the  overall  flux  discrepancy  changes  as  the  partitioning  is  made  finer,  as  represented  by  the 
number  of  sources.  We  found  that  the  use  of  point  sources  makes  a  relatively  small  difference 
(5-7%)  in  the  field  connectivity,  and  that  the  greatest  differences  were  in  connections  which 
contain  relatively  little  flux,  and  thus  are  likely  to  be  unimportant  in  the  energetics  of  the 
region.  Further,  we  found  that  changing  the  side  boundary  condition  typically  results  in  a 
larger  (5-15%)  difference  in  the  connectivity.  Thus  the  assumption  made  by  the  MCC  model 
appears  valid  to  within  approximately  6%,  and  is  not  the  dominant  source  of  variations. 
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Figure  7:  Left:  The  Huxes  of  identical  connections  from  two  different  fields,  MB  (vertical  axis)  and 
PD  (horizontal  axis)  plotted  against  one  another.  Fluxes  are  normalized  to  the  total  ^>tot.  Vertical 
and  horizontal  dotted  lines  show  the  minimum  flux,  0o-  Singlet  connections  (absent  from  one  of 
the  fields)  appear  as  diamonds  below  or  to  the  left  of  these  limits.  For  example,  the  vertical  line 
of  diamonds  are  absent  from  PB  but  present  in  MB  with  fluxes  indicated  by  their  position.  The 
diagonal  dotted  lines  mark  (^PP^^)  (lower).  Statistical  errors 

are  indicated  on  a  few  representative  points.  Plotted  along  the  right  and  top  are  histograms  of  In  ipab 
for  that  field.  The  number  appearing  in  the  histogram  gives  the  total  number  of  connections.  Right: 
The  total  connectivity  difference,  for  different  fields  and  different  partitions.  Comparisons 

are  between  fields  with  different  anchoring  {MB  —  PB)  x,  different  boundaries  {PH  ~  PD)  -f,  and 
differing  in  both  {MB  —  PH)  asterisks.  All  differences  are  plotted  as  a  percentage  of  ^tot-  The 
dashed  lines  shows  {^^\jb~pb  f^^  comparison.  If  the  changes  in  the  connectivity 

from  the  anchoring  and  the  boundaries  were  independent,  one  would  expect  the  dashed  line  to  match 
the  asterisks,  which  is  a  reasonable  approximation  for  most  of  the  range  shown.  [Adapted  from 
Longcope,  Barnes  &  Beveridge  (2009).] 


Comparison  of  Flare  Forecasting  Parameters 

III  recent  years,  the  number  of  parameters  which  have  been  proposed  as  potentially  useful 
for  flare  forecasting  has  proliferated.  One  recently  proposed  parameter  characterizes  the  '‘ef¬ 
fective  connected  magnetic  field”,  Spa,  and  has  been  described  as  “an  efficient  flare-forecasting 
criterion”  (Georgoulis  &  Rust,  2007).  To  test  this  claim,  we  computed  not  just  this  parameter, 
but  several  other  flare  forecasting  parameters  which  have  previously  been  considered,  for  the 
same  data  base  of  active  region  magnetic  field  observations.  The  performance  of  each  of  these 
parameters  was  evaluated  using  a  discriminant  analysis  approach  developed  for  flare  forecast¬ 
ing  (Leka  Sz  Barnes,  2003;  Barnes  et  al.  2007),  as  well  as  the  all-clear  forecasting  approach 
described  in  Georgoulis  Sz  Rust  (2007). 

The  following  parameters  were  use  for  the  comparison: 

=  f  \Bz\  -  the  total  unsigned  magnetic  flux,  often  viewed  as  a  standard  for  judging  other 
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Table  1:  Success  rates  and  skill  scores  for  the  sample  parameters. 


parameter 

success  rate 

skill  score  (binary) 

skill  score  (probability) 

climatology 

0.908 

0.000 

0.000 

^tot 

0.922 

0.153 

0.197 

Ee 

0.916 

0.081 

0.231 

R 

0.922 

0.144 

0.242 

Beff 

0.913 

0.072 

0.220 

parameters. 

Ee  =  J{B^  —  B)^/87r  -  a  proxy  for  the  magnetic  free  energy  (Leka  &:  Barnes,  2003). 

R  =  JfjQpi  \Bi\  -  the  total  unsigned  line-of-sight  field  close  to  high-gradient  polarity-separation 
lines  (Schrijver,  2007). 

~  total  of  the  flux  connecting  sources  divided  by  the  square  of  the 

distance  between  sources  in  an  MCT  model  (Georgoulis  &:  Rust,  2007). 

The  probability  density  of  each  of  these  parameters  for  a  collection  of  over  1200  inagiietograms 
is  shown  in  Figure  8.  In  each  case,  there  is  a  tendency  for  large  flares  to  originate  in  regions 
with  a  large  parameter  \^lue,  but  there  is  a  great  deal  of  overlap  between  regions  which  produce 
flares  within  24  hr  and  those  which  do  not. 

To  quantify  this  behavior,  three  measures  of  success  are  used,  as  shown  in  Table  1.  For 
discriminant  analysis,  a  binary  prediction  is  made  for  each  point,  and  the  success  rate  -  the 
fraction  of  regions  classified  correctly  -  is  used.  Although  the  success  rates  look  impressively 
high,  it  is  important  to  keep  in  mind  that  one  can  obtain  a  success  rate  of  over  90%  by 
predicting  that  all  active  regions  will  remain  flare-quiet  (first  row  of  table).  The  improvement 
over  this  success  rate  for  all  of  the  parameters  is  modest  at  best. 

The  discriminant  analysis  approach  has  also  been  modified  to  produce  probability  forec'asts, 
and  the  climatological  skill  score  is  used  to  judge  their  performance.  The  skill  score  measures 
the  improvement  over  a  uniform  probability  forecast,  with  a  skill  score  of  1  for  perfect  forecasts, 
a  skill  score  of  0  for  no  improvement  over  always  predicting  the  same  probability,  and  a  negative 
score  for  worse  performance.  Again,  the  skill  scores  are  not  dramatically  higher  than  zero, 
consistent  with  the  modest  improvements  seen  in  the  success  rate. 

We  found  that  the  parameter,  proposed  by  Georgoulis  &  Rust  (2007)  does  not  perform 
significantly  better  than  a  similar  parameter,  0tot  considered  by  the  PI  and  Co-I  in  a  previous 
investigation  supported  by  AFOSR.  Indeed,  when  using  discriminant  analysis  to  estimate  the 
success  rate  in  forecasting  the  occurrence  of  major  flares,  it  performs  no  better  than  param¬ 
eters  which  characterize  the  photospheric  magnetic  field,  but  requires  the  additional  step  of 
computing  the  coronal  magnetic  topology.  The  coronal  parameters  may  perform  better  for 
making  “all-clear’  forecasts,  by  predicting  that  regions  with  a  parameter  value  smaller  than 
the  smallest  observed  to  produce  a  flare  will  remain  flare-quiet.  This  approach  is  very  sensitive 
to  a  single  measurement,  and  still  only  successful  in  correctly  predicting  an  all-clear  forecast 
for  a  small  fraction  of  the  possible  all-clear  intervals.  Thus  we  feel  that  pursuing  the  MCC 
approach  is  likely  to  result  in  better  forecasting  results. 

New  and  Refined  Algorithms. 
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Figure  8:  Nonparanietric  density  estimates  for  the  flaring  (green)  and  noiiflaring  (black)  populations 
for  4>tot  (a),  Ee  (b),  R  (c),  and  (d).  The  discriminant  boundary  (50%  probability  forecast)  is 
shown  as  a  vertical  blue  line,  and  the  sample  means  are  shown  as  black/green  vertical  dashed  lines. 
All  the  parameters  exhibit  similar  behavior,  with  a  tendency  for  regions  with  large  parameter  values 
to  be  more  likely  to  produce  an  event,  but  no  clear  separation  between  the  populations.  [Adapted 
from  Barnes  &  Leka  (2008).] 


Magiietogram  Partitioning 

A  critical  step  in  applying  the  MCC  model  is  partitioning  the  active  region  field  into  distinct 
flux  concentrations.  To  be  useful  in  characterizing  active  region  evolution,  it  is  necessary  to 
track  the  evolution  of  these  flux  concentrations  over  a  long  time  interval,  preferably  of  order 
days.  The  approach  we  have  developed  builds  on  a  gradient  based  tessellation  scheme:  one 
partition  consists  of  all  pixels  which  are  strictly  downhill  from  a  local  maximum. 

This  scheme  works  well  on  any  single  magiietogram,  but  when  applied  to  a  time  series,  it 
can  lead  to  large  variations  in  the  partitioning  from  one  time  step  to  the  next.  In  particular, 
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it  is  frequently  the  case  that  a  partition  assigned  to  a  single  maximum  at  one  time  will,  at 
the  next  time,  contain  two  local  maxima  splitting  the  partition  into  two.  To  mitigate  effects 
such  as  this,  we  use  an  approach,  similar  to  Barnes,  Loiigcope  &  Leka  (2005),  of  adjusting 
each  partition  to  minimize  its  differences  with  a  reference.  Barnes,  Longcope  &  Leka  (2005) 
constructed  a  reference  by  averaging  an  approximately  hour-long  time  series  of  magnetogranis 
with  a  cadence  of  a  few  minutes.  For  the  much  longer  time  series  that  may  be  needed  to  apply 
the  MCC  model,  that  approach  will  no  longer  work  because  of  the  significant  evolution  of  the 
photospheric  field. 

Instead  of  a  time-averaged  reference  for  the  entire  time  series,  our  new  method  constructs  a 
new  reference  for  each  magnetogram  by  advecting,  with  the  Local  Correlation  Tracking  (LCT) 
velocity  field,  the  partition  from  a  neighboring  time  step.  For  the  example  presented  here  (see 
Figure  9),  which  contains  regions  of  rapid  flux  emergence,  it  was  preferable  to  start  with  the 
last  time  step,  and  advect  the  partitions  backwards  in  time.  Instances  where  a  new  source 
region  appears  will  seem,  in  the  reversed  time,  to  represent  the  disappearance  of  an  existing, 
labeled  partition. 

For  the  initial  time  the  field  is  first  smoothed  by  performing  a  potential-field  extrapolation 
to  a  height  of  1  pixel.  This  reduces  the  number  of  local  maxima  in  \Bios\  resulting  from  noise. 
The  gradient  based  tessellation  scheme  is  then  applied  to  assign  a  label  to  each  pixel  with  a 
field  strength  exceeding  50  G. 

This  results  in  such  a  large  number  of  partitions  (typically  several  hundred)  in  regions  of 
plage  as  to  be  unwieldy  for  subsequent  calculations.  To  simplify  the  plage  while  maintaining 
the  structure  in  sunspots,  a  saddle  point  merging  is  applied.  That  is,  the  field  is  evaluated  at 
the  saddle  point  between  adjacent  partitions,  and  if  it  is  within  300  G  of  either  maxima  of  the 
two  partitions,  the  smaller  flux  partition  is  assigned  the  same  label  as  the  larger  flux  partition. 
With  this  criterion,  the  plage  is  simplified  to  a  tractable  number  of  sources  (generally  less  than 
100).  Our  hope  in  this  approach  is  that  the  important  effects  in  the  corona  for  solar  flares  and 
CMEs  are  rooted  in  regions  of  strong  field,  so  the  exact  representation  of  weak  fields  is  not 
crucial. 

In  order  to  validate  the  partitioning,  a  comparison  has  been  made  between  the  centroid 
velocity  and  the  flux- weighted  average  of  the  LCT  velocity  for  each  partition.  If  the  partitioning 
algorithm  is  consistently  tracking  the  same  photospheric  flux,  these  velocities  should  be  the 
same.  Figure  9  shows  the  results  of  this  comparison  for  NO  A  A  AR9114.  This  region  was 
chosen  because  the  velocity  of  its  central  sunspot  was  also  the  subject  of  analysis  by  Brown  et 
al.  (2001),  so  that  additional  comparisons  could  be  made.  Both  qualitatively  and  quantitatively, 
there  is  good  agreement  between  the  two  methods  for  determining  the  velocity.  Thus  we  are 
confident  that  the  partitioning  is  indeed  following  concentrations  of  photospheric  flux  as  an 
active  region  evolves. 

Null  Finding  Procedure 

The  null  finding  algorithm  begins  by  sorting  sources  into  ascending  (descending)  order 
for  regions  with  a  net  negative  (positive)  flux.  Neglecting  all  other  sources,  the  null  point 
associated  with  the  first  two  sources  is  determined  analytically.  Each  additional  source  is  then 
reintroduced,  initially  with  a  flux  10“^  smaller  than  its  actual  value.  The  new'  null  point  thus 
introduced  is  located  from  an  initial  guess  given  by  considering  the  new  source  plus  a  uniform 
field  equal  to  the  field  due  to  all  the  other  sources  at  the  location  of  the  new  source.  This  initial 
guess  for  the  null  location  is  refined  using  a  globally  convergent  multi-dimensional  Newlon- 
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Figure  9:  Comparison  of  the  centroid  and  LCT  velocities  for  NOA A  AR  91 14.  Left:  A  high  resolution 
MDI  magnetogram,  with  axes  labeled  in  arc-seconds  from  disk  center.  The  dotted  line  bounds  the 
region  in  which  the  LCT  velocity  is  determined.  For  partitions  entirely  within  this  region,  the  velocities 
as  determined  by  both  LCT  and  fitting  the  location  of  the  centroid  are  shown  as  arrows,  with  length 
proportional  to  the  speed.  There  is  qualitatively  good  agreement  in  the  magnitude  and  direction  of 
the  velocity.  Right:  The  centroid  velocity,  dxajdi^  vs  the  flux- weighted  average  of  the  LCT  velocity 
at  each  pixel,  Ua-  The  x-component  of  the  each  velocity  is  shown  wnth  a  4-;  the  ^-component  is  shown 
with  a  X.  The  dashed  line  is  the  regression  line.  The  correlation  coefflcient  is  p  =  0.77,  indicating 
good  agreement  between  the  methods.  [Adapted  from  Longcope,  Ravindra  &  Barnes  (2007).] 


Rap^hson  root  finding  algorithm.  The  charge  on  the  new  source  is  then  slowly  ramped  up  to 
its  final  value,  and  the  location  of  each  null  point  is  updated  by  using  its  previous  location  as 
an  initial  guess  for  the  Newtoii-Raphson  algorithm. 

At  each  increment  in  the  magnitude  of  the  source  charge,  checks  are  made  for  local  separator 
bifurcations  and  for  pitchfork  bifurcations  both  of  which  result  in  the  creation  of  two  new  null 
points.  This  procedure  is  likely  to  fail  when  more  than  one  bifurcation  occurs  during  a  single 
incrementing  of  the  source  charge.  To  guard  against  this  behavior,  while  also  maintaining  a 
reasonable  speed  for  the  algorithm,  once  each  source  has  been  fully  reintroduced,  the  Euler 
characteristics  given  in  Longcope  &  Klapper  (2002)  are  checked.  If  they  are  not  satisfied, 
then  the  source  is  turned  off  and  reintroduced  more  slowly  until  the  Euler  characteristics  are 
satisfied. 

Technical  Feasibility: 

At  this  stage,  it  is  difficult  to  judge  the  technical  feasibility  of  this  approach  to  predicting 
flare  properties.  We  have  determined  that  applying  the  MCC  model  to  short  time  series  of 
vector  magnetograms  does  not  recover  the  energy  build  up,  but  this  is  almost  certainly  because 
the  process  of  building  up  energy  occurs  on  a  time  scale  of  days.  We  have  developed  most 
of  the  necessary  tools  to  apply  the  MCC  model  to  longer  time  series,  but  an  appropriate 
source  of  data  is  not  presently  available,  although  we  anticipate  that  the  HMI  instrument  on 
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SDO  will  provide  such  data  in  the  near  future,  with  a  launch  data  presently  schedule  for  2009 
November.  Until  these  data  are  analyzed,  it  is  impossible  to  judge  how  well  the  model  will 
work.  Nevertheless,  the  code  which  has  thus  far  been  developed  is  capable  of  running  in  real 
time,  making  it  a  potentially  useful  prediction  tool. 
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